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ABSTRACT 

During the early stages of galaxy evolution, the metallicity is generally low and nearby 
metal-poor star-forming galaxies may provide templates for primordial star formation. 
In particular, the dust content of such objects is of great importance since early molec- 
ular formation can take place on grains. To gain insight into primeval galaxies at high 
redshift, we examine the dust content of the nearby extremely low- metallicity galaxy 
SBS 0335-052 which hosts a very young starburst (< 10^ yr). In young galaxies, the 
dust formation rate in Type II supernovae governs the amount of dust, and by in- 
corporating recent results on dust production in Type II supernovae we model the 
evolution of dust content. If the star-forming region is compact (< 100 pc), as sug- 
gested by observations of SBS 0335-052, our models consistently explain the quantity 
of dust, far-infrared luminosity, and dust temperature in this low-metallicity object. 
We also discuss the H2 abundance. The compactness of the region is important to H2 
formation, because the optical depth of dust for UV photons becomes large and H2 
dissociation is suppressed. We finally focus on implications for damped Lya systems. 
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1 INTRODUCTION 

How dust forms and evolves in primordial galaxies needs 
to be considered before we can understand the chemical 
and thermodynamical evolution of a metal-poor interstellar 
medium (ISM). Dust grains absorb stellar light and reemit it 
in far infrared (FIR)n. Indeed, the FIR spectral range rep- 
resents a unique opportunity to study the dust properties 
and distribution, and dedicated space missions are planned 
in these wavelength bands (ASTRO-Ff|, SIRTE^, HerscheQ 
etc.). Dust grains also drastically accelerate the formation 
rate of molecular hydrogen (H2), expected to be the most 
abundant molecule in the ISM. Hydrogen molecules emit 
vibrational-rotational lines, thus cooling the gas. The cool- 
ing rate of a galaxy consequently depends strongly on the 
dust content and its effect on H2 abundance. 

The production of dust is expected at the final stages 



^ For SBS 0335-052, large energy is emitted in mid-infrared be- 
cause the dust temperature is high (Dale et al. 2001). But in this 
Letter, we use the term "FIR emission" for the thermal emission 
from dust grains. 

^ http://www.ir.isas.ac.jp/ASTRO-F/indcx-e.html 

^ http://sirtf.caltech.edu/ 

* http:/ /astro. estec.esa.nl/First/ 



of stellar evolution, and Type II supernovae (SNe II) are the 
dominant source for the production of dust grains in young 
star-forming galaxies. The formation in stellar winds from 
evolved low-mass stars can also contribute considerably, but 
the cosmic time is not long enough for such stars to evolve 
at high redshift [z). However, dust is also destroyed by SN 
shocks (Jones, Tielens, & HoUenbach 1996). The detailed 
modelling of dust evolution in star-forming regions therefore 
requires an accurate treatment of both types of processes. 

Here we have developed a model for the evolution of 
dust content in primordial galaxies. To compare our model 
with observations, we need age information and the FIR 
energy distribution. For high-z primordial galaxies, present 
experiments cannot determine these quantities accurately 
enough. Nevertheless, because of their low metallicity and 
active star formation, blue compact dwarf galaxies (BCDs) 
may be viable candidates for nearby primeval systems. We 
have focused on the evolution of dust content in one of these, 
SBS 0335-052, since this object may be experiencing its first 
burst of star formation (age < lO'' yr; Vanzi et al. 2000). 

In this Letter, we first model the dust content of a young 
galaxy whose age is less than 10* yr (§ |2[). Then, we com- 
pare the model predictions with the observed quantities of 
the extremely low-metallicity BCD, SBS 0335-052 (§ ||). Fi- 
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nally, we comment on implications for damped Lya systems 
(DLAs). 
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2 MODEL DESCRIPTION 

2.1 Evolution of Dust Mass in a Young Galaxy 

We aim at modelling the evolution of dust content in a young 
galaxy whose age is less than 10* yr. We focus on the com- 
parison with the extremely metal-poor (1/40 Zq) BCD SBS 
0335-052, one of the few possible local counterparts of a pri- 
mordial galaxy in the nearby Universe (Vanzi et al. 2000). 
Because of its young age, the contribution of winds from 
late-type stars to dust formation is assumed to be negligible. 
We consider the dominant (brightest) star-forming complex, 
and assume that grains produced by SNe II are accumulated 
within this region. This assumption is consistent with high- 
resolution infrared imaging observations which reveal that 
the distribution of dust is as compact as the size of the 
star-forming region (Dale et al. 2001; Hunt, Vanzi, & Thuan 
2001). 

We start by considering a generic star-forming region 
in a galaxy. Our model is applicable to young (< 10* yr) 
bursts. Because of its young age, the contribution of Type 
la supernovae to the total supernova rate is neglected. The 
rate of SNe II is given by 



lit) 



(m) dm , 



(1) 



where 7/)(f) is the star formation rate (SFR) at t (we de- 
fine t = at the beginning of the star formation), 0(m) is 
the initial mass function (IMF), Tm is the lifetime of a star 
whose mass is m, and we assumed that stars with m > 8 A^q 
produce SNe II. In this Letter, we assume a constant SFR, 
Tp ~ ipo, and a Salpeter IMF {ip{rn) oc m~^ "^^; the mass 
range of stars is 0.1-60 M0) to obtain a first estimate. Our 
assumption of constant star formation is applicable if the 
time from the starburst is small. Indeed, the difference be- 
tween the instantaneous burst and the constant star forma- 
tion history is small on small timescales (< 10^ yr, the age 
of SBS 0335-052; Vanzi et al. 2000). 

In general, dust destruction by SNe can be important, 
but for SBS 0335-052 it is negligible for the following reason. 
The destruction timescale tsn is estimated to be (McKee 
1989; Lisenfeld & Ferrara 1998) 



TSN 



7eMs(100 km s"^) 



(2) 



where Mg is the total gas mass in a star-forming region, 
Ms(100 km s"^) = 6.8 x 10^ Mq (Lisenfeld & Ferrara 1998) 
is the mass accelerated to 100 km s~^ by a SN blast, 7 
is the SN II rate, e ~ 0.1 (McKee 1989) is the efficiency 
of dust destruction in a medium shocked by a SN. We as- 
sume the relation between stellar mass and lifetime as shown 
in Inoue, Hirashita, & Kamaya (2000a). The SN II rate 
increases with time, and reaches a constant value, if the 
star formation rate is constant in time. Numerically, we find 
that 7 (yr^^) ~ 3 x 10"^ Vo [Mq yr"^) around the age of 
5 Myr. We also assume that the star formation efficiency 
is the order of 5% as known empirically from observations 
(e.g., Inoue, Hirashita, & Kamaya 2000b). This means that 
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Figure 1. Time evolution of dust amount in a galaxy {solid line). 
The initial star formation rate is assumed to be 1 Mq yr~ ^ . Tlie 
horizontal dotted lines show the value determined observationally 
by Dale ct al. (2001) and the upper limit by Hunt ct al. (2001). 



Mg/ipo ~ 20tmax, where tmax is the possible duration of the 
star- forming activity. Thus, we finally find that tsn 3> tmax- 
It is probable that t < tmax during the ongoing starburst, 
and we can safely conclude that the dust destruction can be 
neglected. 

Therefore, the only contribution to total dust mass in 
the galaxy (Md) comes from the supply from SNe II. Then 
the rate of increase of Md is written as 



Md = rudj , 



(3) 



where md is the typical dust mass produced in a SN II. To- 
dini & Ferrara (2001) showed that rrid varies with progeni- 
tor mass and metallicity. There is also some uncertainty in 
the explosion energy of SNe II. The Salpeter IMF-weighted 
mean of dust mass produced per SN II for the 1) Z = 0, 
Case A,2)Z = 0, Case B, 3) Z = W'^Zq, Case A, and 4) 
Z = W'^Zq, Case B are 1) 0.22 Mq, 2) 0.46 Mq, 3) 0.45 
Mq, and 4) 0.63 Mq, respectively {Z is the metallicity, and 
Cases A and B correspond to low and high explosion en- 
ergy|^, respectively). We adopt the average of the four cases, 
i.e., rud — 0.4 Mq, but the calculated dust mass is propor- 
tional to TTid. Therefore, there is uncertainty in md and Afd 
by a factor of ~ 2. Since we are considering the first burst 
of star formation, we assume Aid = at t = 0. Because dust 
production rate is proportional to 7 which is proportional 
to the star formation rate ipo, Md calculated by our model 
directly scales as tpo- 

In Figure |l|, we show the dust mass evolution for a star 
formation rate of 1 Mq yr~^, a typical value for SBS 0335- 
052 (Hunt et al. 2001). In Figure h, we also show the values 
of dust masses derived by Dale et al. (2001) and Hunt et al. 
(2001). The dust mass in Hunt et al. (2001) was calculated 
with a unit filling factor, and is considered to be an upper 
limit. 
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2 X 10^^ ergs for Case A and Case B, respectively. 
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Figure 2. Time evolution of tfie far-infrared (FIR) luminosity 
for the same star formation history as in Fig. hi. The radius of 
the galaxy is assumed to be 300, 100, and 30 pc (solid, dotted, 
and dashed lines, respectively) The observed FIR luminosity is 
shown by the horizontal dot-dashed line (~ 5 X 10* Lq; Dale et 
al. 2001). 



2.2 Dust Temperature and FIR Luminosity 

Probably the most direct way to measure dust content is to 
observe its FIR emission. We now derive the evolution of 
FIR luminosity and dust temperature. Because of the large 
cross section of dust against ultraviolet (UV) light and the 
intense UV radiation field in a star-forming galaxy, we can 
assume that the FIR luminosity is equal to the absorbed 
energy of OB stellar radiation (~ UV radiation). First, we 
estimate the optical depth for the UV radiation as follows: 

_ 2^ 9 QuvMd 

T-dust Tra Qvyn^rsF = — — , (4) 

IbTT adrgp 

where a is the grain radius (spherical grains are assumed), 
na^Quy is the typical absorption cross section for UV light, 
rid is the number density of grains, rsF is the radius of 
the star-forming region (a spherical star-forming region is 
assumed), and S is the grain material density. We adopt 
a ~ 0.03 Atm (Todini & Ferrara 2001), Quv ~ 1, and 5 = 2 
g cm~^ (Draine & Lee 1984). Then we obtain the energy 
absorbed by dust. We assume that all the absorbed energy 
is reemitted in the FIR. Thus, the FIR luminosity Lfir be- 
comes 



Lfib. = Lob[1 - exp(-rdus 



(5) 



Lob is calculated assuming a mass-luminosity relation of 
OB stars (> 3 Mq; we adopt the fit by Inoue et al. 2000a). 

In Figure ^ we show the evolution of FIR luminosity 
for rsF = 300, 100, and 30 pc (solid, dotted, and dashed 
lines, respectively). Dale et al. (2001) have estimated that 
rsF ^ 80 pc for SBS 0335-052. We reexamine this finding to 
check the consistency with the dust supply rate from SNe 
II. The observed FIR luminosity (~ 5 x 10* Lq) for SBS 
0335-052 is reached in 8, 3, and 1.2 Myr for rsF = 300, 100, 
and 30 pc, respectively. 

The dust temperature Td is determined from the emis- 
sivity given by Draine & Lee (1984); i.e., Qabs(A) = 3 x 



Figure 3. Time evolution of the dust temperature for the same 
star formation history in Fig. |l[ The radius of the star-forming 
region is assumed to be 300, 100, and 30 pc (solid, dotted, and 
dashed lines, respectively). The observed dust temperature is 
shown by the horizontal dot-dashed line (80 K; Dale et al. 2001). 



10^*(A/100 fim)^^, where (5abs(A) is the absorption effi- 
ciency (the definition is the same as Draine &i Lee) as a 
function of wavelength A. Then, we finally obtain the fol- 
lowing relation: 



Td = 20 



Lfih/Lq 



2.5 X lO^A'/d/Me 



1/6 



K. 



(6) 



In Figure ^, we show the evolution of dust temperature. 
The predicted temperature (~ 45-120 K) is higher than that 
observed in normal spiral galaxies (~ 20 K), but such a high 
temperature (~ 80 K; dot-dashed line; Dale et al. 2001) 
is observed in SBS 0335-052. However, the observed FIR 
luminosity and temperature are uncertain because the peak 
of FIR spectrum has not been detected yet. Thus, the three 
quantities (the size, FIR luminosity, and dust temperature) 
are simultaneously accounted for by our model. Such a high- 
dust temperature is also predicted for high-z objects (e.g., 
Totani & Takeuchi 2001). 

2.3 H2 Abundance 

The abundance of H2 affects the cooling rate of gas. In par- 
ticular, to obtain insight into the molecular cooling of a pri- 
mordial system, nearby star-forming dwarf galaxies consti- 
tute an important "laboratory". Thus, we next calculate the 
evolution of molecular hydrogen fraction. 

We assume that molecules form on dust grains and are 
destroyed (dissociated) by the UV field (both processes oc- 
cur on a much shorter timescale than the lifetime of OB 
stars; Kamaya & Hirashita 2001). Thus, we evaluate the 
two rates and assume equilibrium. We adopt the formation 
rate of H2 by HoUenbach & McKee (1979): 



R = 0.5nH(l - /H2)nd7ra ^^(r) , 



(7) 



where tih is the number density of hydrogen nuclei, is 
the number fraction in molecules (i.e., /h2Wh/2 is the num- 
ber density of H2), v is the mean thermal speed of hydrogen. 
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and S{T) is the sticking coefficient of hydrogen atoms. The 
latter is given by (Omukai 2000) 

^(r) = [1 + o.o4(r + rd)"-'"' + 2 X io-^r + 8 X io~®T^]"^ 

X [l + exp(7.5x 10'(l/75-l/rd))]-\ (8) 

where T is the gas temperature. We take a — 0.03 fim (To- 
dini & Ferrara 2001) and T = 100 K. The gas temperature 
is at present poorly constrained both by absorption line and 
21-cm H I observations; however it should represent a rea- 
sonable estimate for the temperature of the cold, H2-forming 
phase of the Galaxy ISM. 

On the other hand, the dissociation rate /tihs is esti- 
mated empirically by 

^"H^ - 2xl0-"f--— ) 

\2 X 10^^ erg s^^ cm^^y 

x/Haf^H cm"^ s"^ , (9) 

where /uv ~ Lqb exp(— rdust)/7'sF is the UV radiation field 
attenuated by dust. We have assumed for normalization that 
the UV radiation intensity and dissociation rate at the solar 
neighbourhood are 2 x 10""^ erg s~^ cm~^ and 2 x 10~^^ s~^ 
(Jura 1974), respectively, and that 7nH2 linearly scales with 
7uv- Equating the two rates, we obtain /hj as a function of 
age of burst. For a given amount of dust, /hj is indepen- 
dent of the hydrogen number density. Thus, we examine the 
dependence on rsp. Since we assume that Md = at t = 
and consider only the reaction on grains, the initial molec- 
ular fraction is assumed to be zero (/hj = at t = 0) . If H2 
formation in the gas phase is allowed, > at t = 0. In 
this sense, our estimate of /hj is a lower limit. 

The result is shown in Figure ^ for rsp = 300, 100, 
and 30 pc (solid, dotted, and dashed lines, respectively). 
The molecular fraction rises because of the shielding of UV 
radiation by grains. Thus, the dust enrichment by SNe II 
causes a significant shielding of UV radiation field. Once 
the column density of H2 becomes larger than 10^* cm~^, 
the self-shielding effect becomes important (e.g., Draine 
& Bertoldi 1996). The value corresponds to /hj — 3 x 
10"^(nH/100 cm-3)-i(rsF/100 pc)-^ The emission line di- 
agnostics by Izotov et al. (1999) and Vanzi et al. (2000) indi- 
cate that nu ~ 500 cm~'^ in the star-forming region of SBS 
0335-052. With nn = 500 cm"^, the self-shielding value of 
/h2 is achieved at the age of 10, 5, and 2 Myr for rsp — 300, 
100, and 30 pc, respectively. 



3 DISCUSSION 

3.1 Physical Conditions of SBS 0335-052 

For the dust mass. Dale et al. (2001) derived the value of 
Md = 2400 Mq from black-body fitting of the observed 
spectrum of SBS 0335-052. Hunt et al. (2001) derived an 
upper limit of 10^ Mq from the extinction. A starburst of 
age 3 X 10^-5 x 10^ yr reproduces those values (Fig. ^). This 
range is consistent with the age observationally determined 
by Vanzi et al. (2000). For the dust temperature, it is nec- 
essary that the radius of the star-forming region be smaller 
than 100 pc to explain the 80 K component revealed by ISO 
(Dale et al. 2001), although the temperature may be overes- 
timated because of the contamination of the emission from 
very small grains heated stochastically. 
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Figure 4. Time evolution of the molecular fraction for the same 
star formation history in Fig. |l[ The radius of the star-forming 
region is assumed to be 300, 100, and 30 pc (solid, dotted, and 
dashed lines, respectively). 

As observational information on the H2 abundance, 
near-infrared emission lines are available. The line inten- 
sity is dependent on the excitation mechanism and thus the 
quantity of H2 is poorly constrained. However, the detection 
of the H2 lines (Vanzi et al. 2000) implies that there is an ef- 
fective mechanism that shields the UV-dissociative photons. 
The compactness as observed by Hunt et al. (2001) and Dale 
et al. (2001) may explain the existence of H2 for this object 
(Fig. M) because of the shielding effects. 

Finally, we conclude that if the star forming region of 
SBS 0335-052 is compact (rsp ^ 100 pc) and the gas density 
is large (nn > 10^ cm~^), the young age, the dust amount, 
the FIR luminosity, dust temperature, and the H2 abun- 
dance are mutually consistent and explained by our models. 



3.2 Implications for High-Redshift Galaxies 

The main conclusion drawn in this Letter is that star for- 
mation activity in extremely metal-poor systems can suffer 
from extinction. Thus, FIR observations of primeval galaxies 
are important to observe such dust-enshrouded star forma- 
tion activity. 

Our results for the H2 abundance are interesting in con- 
nection with the recent observations of DLAs. There are 
several DLAs whose abundances of dust and H2 are obser- 
vationally derived. Levshakov et al. (2000) assumed equilib- 
rium between the photodissociation of H2 by the background 
UV radiation and formation on grains. They suggested that 
the small dust number density (as low as 10"'^ times the 
Galactic value) can account for the observed small molecu- 
lar fraction (/hj ~ 4x 10~*). We have shown that such a low 
molecular fraction can be due to the dissociation caused by 
recently formed OB stars. We also have suggested that if the 
region is compact and dense, the dust produced by SNe II 
can shield the UV radiation. This can produce a strong cor- 
relation between dust amount and H2 fraction. Indeed, a cor- 
relation between the two quantities has been observationally 
verified. Some DLAs whose depletion is large show molec- 
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ular fractions up to 0.2 (Ge, Bechtold, & Kulkarni 2001). 
Such high values are possible when dust shielding of the UV 
radiation field is effective. 
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